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SUMMARY 

This paper presents the results of a study to determine the accuracy of 
the three dimensional sdid elements available in NASTRAN for 
predicting displacements. Of particular interest in the study is 
determining how to effectively use solid elements in analyzing thick 
optical mirrors, as might exist in a large t^scope. Surface 
deformations due to thermal and gravity loading can be significant 
contributors to the determination of the overall optical guality of a 
telescope. The study investigates most of the soUd elements currently 
available in either COSMIC or MSC NASTRAN. Error bounds as a 
function of mesh refinement and element aspect ratios are addressed. 
It is shown that the MSC solid elements are, in general, more accurate 
than their COSMIC NASTRAN counterparts due to the specialized 
numerical integration used. In addition, the MSC elements appear to be 
more economical to use on the DEC VAX 11/780 computer. 

INTRODUCTION 

Optical mirrors for spaceborne telescopes often have thicknesses which 
are as much as 10 percent of their diameter. Because of this, and also 
due to the need to obtain accurate predictions of optical surface 
deformations in the presence of temperature gradients that are often 
nonlinear, these mirrors are modeled with three dimensional, sdid, 
finite elements. In addition to thermal deformation being a significant 
design driver for these mirrors, deformation due to "gravity release" 
must also be considered. When fabricated and tested, the mirror is in 
a 1— g environment which is "released" once the telescope reaches the 
0-g environment in orbit. 

A common design practice is to fabricate these mirrors from low 
coefficient of thermal expansion (GTE) materials but to mount them to 
conventional metal housings with much larger CTE's. In order to avoid 
stressing the mirror when it is subject to simple bulk temperature 
changes, the mirror is often mounted to its housing through the use of 
(nearly) kinematic mounts that may connect to the mirror in only three 
locations. 

The purpose of this study is to determine how to effectively model such 
mirrors for thermal and gravity loading using the av ai l a ble solid 
elements in NASTRAN. In a study such as this there are three 
significant questions to consider. 

1. How fine a mesh must be used to obtain displacements of the 
optical surface within a certai n error range? 

2. Do the individual scQid finite elements exhibit sensitivity to aspect 
ratios significantly different than 1.0? 
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3 . 


Which of the avai^ble Aments are the most economical to use 
(accuracy vs cost in computer resources) ? 


Questions 1 and 2 above can best be answered in an investigation that 
compietely separates the two effects. That is, the mesh refinement 
smoies should be done using elements with an aspect ratio of 1.0. 
Then, once a fine enough mesh has been reached such that the errors 
are small, the effects of aspect ratio can be investigated by keeping the 
mesh size constant and varying the overall dimensions of the problem, 
thus resulting in each element aspect ratio changing. Obviously, in 
order to accomplish this latter step there must be a theoretical (or some 
^er good comparison) solution to the problem with which to compare 
^e finite element model results since at each step a problem of 
different dimensions (and therefore different theoretical solution) is 
eing modeled. This technique has been used in previous studies for 
two dimensional membrane elements, [1] , [2] . 


With the above conaderations in mind, the sample problem used in the 
Study ^ (described in more detail below) is a cubic slab of ecrual 
dimension in the X-Y plane and whose thickness varies between 
one-twentieth and one-halE of the X-Y plane dimensions (see Figure 1) 
The sMbs were constrained in a kinematic fashion simulating the type of 
repaint often used in actual practice. For the thick slab, a mesh 
refinement can be made of elements having aspect ratios of 1.0 without 
requiring an exorbitant number of elements in the X-Y plane when a 
reasonable number is used through the thickness. For a linear 

through the slab there exist theoretical solutions. 
Thus the^^k ^b, although it does not look much Hke a mirror, is a 
good can^date for the mesh refinement studies when the temperature 
gra^nt loading is used. Also, since the theoretical solution exists for 
arbitrary thickness of the slab, the sample problem is also suitable for 
an ^ment aspect ratio study. Keeping the mesh refinement constant 
1 number of elements in the model) , simply by varying the 

slab thickness, the element aspect ratios must change. 


^ gravity loading the situation is not as good. There is no 
theoretical solution to compare to so it is not practical to attempt the 
aspect rato study with this loading. However, if the various elements 
used in the study show a trend toward convergence to an answer as 

• ^ ^ refined, the mesh refinement studies can provide useful 

information with this loading. useiui 

PROBLEM DETERMINATION 

^ shows the geometry, coordinate system, boundary conditions, 
vfrim ^ information used in this study. The constraints are 

J^ema-^ and the problem is symmetric about the x=0 plane. That is 

S!?, \ ^ plane. Using this constraint,' 

only half the slab was modeled for the study. 

elements were tested in this study. This 
ind-uded hexahedr^ elements with up to two midside nodes and wedge 

shaped ^ments with up to one midside node. Figure 2 displays the 
different types examined. cue 
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The mesh geometry is shown in Figure 3. AH mod^ had half the 
number of elements through the thickness as in the other two 
dimensions. Only the diamond pattern shown in Figure 3 was used for 
wedge shaped elements and the triangular shape was always in the x-y 
plane as recommended in [4] . When an element aspect ratio was 
desired, the value of t was changed while the number of elements 
through the thickness remained unchanged. This created shorter or 
squat elements. 

Temperature Gradient Loading 

The temperatures applied to the test model varied only in the z 
direction. A linear gradient was created which ranged from +h to -h 
degrees centigrade. The equation for this gradient is: 

T^'/zO-Vt) 


T* tempere»^or<. , thickness 

An exact theoretical solution exists for a linear temperature gradient 
[3] and in this case the equations are: 

u - U -t Uo V ^ 

V == V -f So - £©k. t X0e. 

3 \a) + Ulo - X©yo + 

where 

U,V, w a d'lsplacemerd^ i = GTE 

V * ^(Yz” S - 

u) “ z ■* Vz-t 0 ^ + -1 

Using the boundary conditions as shown in Figure 1 we get: 

u = 

v = «<.<Y2 .- Vt)y 

u; •= ot\(Va-M,'ia -v Vs.-tC'i?' ■v'1% 

Gravity Loading 

Since many telescopes operate in a zero gravity environment, surface 
deformations due to their original one G environment are important in 
evaluating their optical quality. The load used in this study was 

28 



applied in the -z direct±)n wibh the supports at the bottom surfece. No 
theoretical solution is available for this case. 


A thin slab (i.e.^ t=.0504M (2")) model was used for the gravity 

loading. Since this created Aments with aspect ratios of 5, only 
elements which showed no susceptibility to aspect ratio errors in the 
temperature load case were used. 

RESULTS 


Temperature Gradient Loading 

For the linear temperature gradient loading the results of the mesh 
refinement and aspect ratios studies are shown in Figures 4 and 5. 
Figure 4 shows the error in displacement at pcdnt A as a function of 
mesh size (N) for each of the nine scQid elements in the study. The 
number of elements in any one model is N x N x N/2 in the x, y, z 
directions, respectively, for the complete problem and are half this 
amount for the one— half of the slab actually modeled due to symmetry. 
In order to keep the element aspect ratio 1.0 for this part of the 
study, the slab dimensions were chosen to be .508 x .508 x .254M (20 
X 20 X 10"). 

The linear temperature gradient, coupled with the stress-free mounting 
used in the problem, produces a Unear strain variation through the 
tliickness of the slab. It is to be expected, therefore, that 
3Soparametric elements of high enough order, when integrated using 
Gauss quadratur with a sufficient number of pcdnts, will yield exact 
results (even when the problem contains only one element through the 
thickness). This is true for the COSMIC quadratic and cubic elements, 
CIHEX2 and CIHEX3. In addition, it should also be true for the MSC 
elements CHEXA (20 node) and CPENTA (15 node) which are also 
quadratic isoparametric elements but which use reduced integration for 
selected terms in the stiffness matrix. Reduced integration is used in 
conjunction w^ many isoparametric elements in an effort to improve on 
the overly stiff behavior of some of the lower order of these elements 
[5] -[8]. When used in situations in which there is primarily bending 
behavior, the lower order isoparametric dements have a "parasitic" 
shear introduced due to the allowable modes of deformation as defined 
by the polynomi^ displacement field . This parasitic shear can be 
removed by sdectLvely reducing the order of numerical integration used 
an generating the element stiffness matrices. The resxilt is an element 
wh^h no longer demonstrates interelement displacement compatibility but 
which, paradojdcalLy , often performs better. Table 1 lists the order of 
the numerical integration used for the isoparametric elements used in 
this study. The exact det ails of the reduced integration for the MSC 
elements is not made dear in their documentation but is believed to be 
such that the shear terms are represented by only one point for the 
linear displacement polynomial elements. 

Figure 4 shows that the MSC linear isoparametric dements give exact 
behavior as wdl as the higher order elements which were expected to 
be exact. The linear elements, hke their COSMIC counterpart CIHEXl, 
are based on products of linear shape functions which are capable of 
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representing all constant states of strain and some linear states as well. 
However, the linear normal strains in the x, y directions are not 
represented by these elements so that the exact behavior of the 8 node 
CHEXA MSC element is not fully understood. It should be anticipated 
however, that the 8 node CHEXA would be better for this temperature 
loading than the COSMIC CIHEXl due to the reduced integration used. 

In order to determine the sensitivity of these elements to element aspect 
ratios different than 1.0, only the finest mesh for each element was run 
with the slab thickness decreasing in successive runs from .254 M (10") 
to .0254 M (1"). This resulted in element aspect ratios changing from 
1.0 to 10 (that is, the element dimensions in the x-y plane, in ratio to 
its thickness changed from 1.0 to 10). Figure 5 shows how the 
displacement errors are affected by these changes in element aspect 
ratio. The only isoparametric element which shows sensitivity to aspect 
ratio is the COSMIC linear element CIHEXl. 

Since the slab thickness has decreased from .254 M (10") to .0254 M 
(1") while the element aspect ratio increased from 1 to 10, the slab 
deflections are due more to bending now than they were when the 
thickness was .254 M (10"). As mentioned previously, the CIHEXl has 
a parasitic shear due to bending so that it might be expected that the 
deflection errors would worsen as the slab exhibits more bending . T he 
other isoparametric elements avoid this problem, to some degree or 
another, by either having higher order displacement polynomials and/or 
reduced integration. 

Gravity Loading 

As with the temperature loading study , the gravity loading study began 
with the thick slab in order to preserve the 1.0 element aspect ratio 
during the mesh refinement portion of the study . As mentioned 
previously, there is no theoretical solution to a thick slab loaded by a 
body force when constrained in the kinematic fashion used in this 
study. Thus, Figure 6, which shows the mesh refinement results, has 
actual displacement (rather than displacement error) plotted versus 
mesh size. The results shown on Figure 6 indicate no convergence has 
taken place for the range of mesh sizes used. Apparently, the reason 
for no convergence is due to the concentrated loads at the support 
points (along either side of the slab) . The mesh for this kind of a 
problem really needs to be more refined in these areas rather than just 
using a uniform mesh for the whole problem. In addition, the "point" 
support must be distributed over a finite surface area to avoid the 
singularity (in displacement and stress) that exists, in the limib, as the 
mesh size (N) becomes infinitely large. However, in a three 

dimensional problem, this mesh refinement is very difficult to do without 
the aid of some automatic mesh generation program (which was not 
available for a this study) . 

The difficulty in convergence evidenced by the results in Figure 6 was 
thought to be compounded by the fact that the slab is q\iLte thick (one 
half of the in-plane dimensions) . Since the overall study is really 
aimed at how to model moderately thick mirrors it was decided to change 
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the slab to one that had a thickness of one-tenth the in-plane 
dimensiDns but keeping the element mesh arrangement the same (N x N 
X N/2). Since this would resvilt in elements with aspect ratios of 5, 
only those elements that exhibited no aspect ratio sensitivity in the 
thermal load study were used. This included the 1 and 2 mid-side node 
COSMIC isoparametric elements and aH of the MSC isoparametric 
elements. 

Figure 7 shows the results of the mesh refinement study using the thin 
dab (.504 X .504 x .0504 M (20 x 80 x 2")). Again it is obvious that 
convergence ^ has not been reached by any of the element types even 
though the finest mesh for the MSC 20 node hexahedral was 10 x 10 x 5 
and contained 4064 degrees of freedom for the half model. This run 
. took nearly six hours of CPU time on the VAX 11/780. An attempt to 
run this element for the 12 x 12 x 6 mesh foiled due to the size of a 
matrix multiply requiring more memory than was available. Table 3 
shows how large these problems can become (for the half model) using 
the N X N X N/2 uniform mesh. The run that failed contained nearly 
6700 degrees of freedom. 

From the lack of convergence for this more practical problem it is 
obvious that a uniform mesh size is very impractical to use for a 
loading and constraint system like those encountered for some large 
mirrors. The local deformation in the vicinity of the support pcdnts 
requires refinement. Figure 8 illustrates this by showing the 
displacement at point B which is directly above one of the support 
points (which carries one half the weight of the slab) . These 
displacements obviously are much smaller than those at point a, 
however, they are still growing significantLy as the mesh is refined. In 
an effort to see if the convergence problem is really associated with the 
concentrated loading at the support pcdnts, the relative deflection 
b^ween points A and B was determined. Figure 9 shows this relative 
displacement and it appears to be converging indicating that the front 
surface (i.e. the z = t surface) shape is defined but that the whole 
slab is still "sinking” over the point support under B (and its 
symmetric counterpart at y = z = 0 and x = 1/2) . Figure 9 is really 
only u^ful in demonstrating that the mesh really needs r ef ining (and 
the point support needs to be distributed over a finite area) in the 
vicinity of the supports in order to obtain useful information about the 
deflections of the slab. 

Timing Study 


Table 2 shows all of the elements tested and the computer time required 
to generate their stiffness matrices on the DEC VAX 11/780. The 
COSMIC isoparametric elements with one and two midside nodes take an 

appreciable amount of time and limit the practical problem sizes which 
can be run. 

In realty, it is the total computer costs for a given accuracy that is of 
most imporfonc^. One of the original goals of this study was to 
address this issue. However, due to the difficulty with obtaining 
convergence for the gravity loading , a "comparison" answer was not 
able to be found with which to guage the absolute accuracy of the 
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gravity loading results. This precluded obtaining a plot of accuracy 
versus cost for this loading. However, based on the solution times for 
the data of Figures 7-9, it appears that the MSC hexahedra with one 
midside node is the best element in terms of accuracy versus cost. 

Conclusions 


Most of the available isoparametric soli d elements in both COSMIC and 
MSC N AST RAN appear to be well suited for thermal deformation 
analyses of kinematically supported mirror-type structures. An 
exception to this is the COSMIC 8 noded hex element (CIHEXl) which 
exhibits large errors when the element aspect ratio deviates from 1.0; 
the element should not be used for aspect ratios above 2.0. 

For gravity loading, the situation is much less dear. Use of a uniform 
mesh in three dimensional problems is an attractive choice due to the 
complexity of envisioning nonuniform meshes in more than two 
dimensions. However, the study indicates that with point supports, an 
extremely fine mesh woiold have to be used to get convergence for any 
of the elements used and the "point" support must be distributed over 
some fixed fini te area as the mesh becomes more and more refined. It 
appears that mesh refinement in the vicinity of the supports is the only 
realistic way to model these structures and use of some automatic mesh 
generation program to accomplish this is a necessity. 

From the standpoint of economy, the MSC 20 node hex element (CHEXA) 
appears to be by far the best of all of the available elements. 
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Table 1 


Gauss Quadrature Numerical 
IntegratiDn Mesh Used For 
Isoparametric Elements 

N AST RAN Isoparametric Number Displacement Gauss QuadratLve 

Version Element of Nodes Polynomial Mesh Points 

COSMIC CIHEXl 8 linear 

CIHEX2 20 quadratic 

CIHEX3 32 cubic 

MSC CHEXA 8 linear 

CHEXA 20 quadratic 

CPENTA 6 linear 

C PENT A 15 quadratic 

(1) AH MSC dements use a numerical integration for selected terms 
(shear) at a number of Gauss points reduced from those indicated. 
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Tab]fi 2 



Element Generation 

EMG + EMA time 

NASTRAN 


per element 

Version 

Element 

(C.P.U. sec) 2 

COSMIC 

CHEXA2 

2.1 


CWEDGE 

1.0 


CIHEXl 

1.2 


CIHEX2 

10-13 


CIHEX3 

23-29 

MSC 

CHEXA (8 nodes) 

.8 


CHEXA (20 nodes) 

4.5 


CPENTA (6 nodes) 

.4 


C PENT A (15 nodes) 

1.8 


* These times were taken from models run on a VAX 11/780 using 
COSMIC versions 17.6 and April 1983 and MSC version 63. 
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Table 3 


Number of Degrees of Freedom for Specdfic Model Sizes 


N 


HEX 

or WEDGE 
0 

midside nodes 


HEX 

1 

midside node 


WEDGE 

1 

midside node 


HEX 

2 

midside nodes 


2 

27 

80 

92 

133 

4 

117 

383 

455 

649 

6 

305 

1052 

1268 

1799 

8 

627 

2231 

2711 

3835 

10 

1119 

4064 

4964 

7009 

12 

1817 

6695 

8207 

11,573 
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APPENDIX - NOTATION 


A^£ = ELEMENT ASPECT RATIO (le/te) 
^ = PROBLEM SLAB LENGTH 

= ELEMENT LENGTH 
E = YOUNGS MODULUS 
t, = PROBLEM SLAB THICKNESS 
“te = ELEMENT THICKNESS 
U;V,U; = DISPLACEMENTS 

= COEFFICIENT OF EXPANSION 
p = DENSITY 
V = POISSON RATIO 
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FIG. 1 

TEST PROBLEM 



♦ INDICATES DEGREE OF FREEDOM CONSTRAINED FOR 
KINEMATIC MOUNTS 

v=w = 0ATY = Z = 0, X = ±«/2 
u=w = 0ATX=Z=0, Y =f/2 

A = (0.,0.,t) 

B = (</2,0.,t) 

MATERIAL INFORMATION (ALUMINUM) 

E = 6.89 X N/M^ (10 x 10® LB/IN^) 
p = 2.71 X 10® KG/M® (.098 LB/IN®) 
a = 22.7 X 10'® /K 
V = .33 
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FIG. 2 


ELEMENT TYPES 



8 NODES 



CHEXA2 10 TETRAHEDRA (COSMIC) 

CIHEX1 LINEAR ISOPARAMETRIC (COSMIC) 
CHEXA LINEAR ISOPARAMETRIC (MSC) 


CIHEX2 QUADRATIC ISOPARAMETRIC (COSMIC) 
CHEXA QUADRATIC ISOPARAMETRIC (MSC) 



CUBIC ISOPARAMETRIC (COSMIC) 



CWEDGE 3 TETRAHEDRA (COSMIC) 
CPENTA LINEAR ISOPARAMETRIC (MSC) 


15 NODES CPENTA QUADRATIC ISOPARAMETRIC (MSC) 
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FIG. 3 


MESH GEOMETRY 




WEDGE PATTERN 
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FIG. 4 


Z DISPLACEMENT ERROR AT (A) 
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FIG. 5 

Z DISPLACEMENT ERROR AT (A) 
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LINEAR TEMPERATURE GRADIENT 
(ASPECT RATIO STUDY) 
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FIG. 6 

Z DISPLACEMENT AT (A) 

THICK SLAB — GRAVITY 
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FIG. 7 

Z DISPLACEMENT AT (A) 


THIN SLAB-GRAVITY LOADING 
(MESH SIZE STUDY) 
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FIG. 8 
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Z DISPLACEMENT AT (B) 

THIN SLAB-GRAVITY LOADING 


(MESH SIZE STUDY) 
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FIG. 9 

RELATIVE DISPLACEMENT (A-B) 

THIN SLAB-GRAVITY LOADING 
(MESH SIZE STUDY) 
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k. CPENTA 6 NODES 
CPENTA 15 NODES 
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MESH SIZE, N = Vce 




